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The polarization analysis of the light is typically carried out using modulation 
J"}" 1, schemes. The light of unknown polarization state is passed through a set of 

known modulation optics and a detector is used to measure the total intensity 
CNj ■ passing the system. The modulation optics is modified several times and, 

i— 1 1 with the aid of such several measurements, the unknown polarization state 

of the light can be inferred. How to find the optimal demodulation process 

q . has been investigated in the past. However, since the modulation matrix 

+-> ■ has to be measured for a given instrument and the optical elements can 

, &, • present problems of repeatability, some uncertainty is present in the elements 

of the modulation matrix and/or covariances between these elements. We 

analyze in detail this issue, presenting analytical formulae for calculating the 
00 \ 
O \ covariance matrix produced by the propagation of such uncertainties on the 

demodulation matrix, on the inferred Stokes parameters and on the efficiency 
■<^j- . of the modulation process. We demonstrate that, even if the covariance 

matrix of the modulation matrix is diagonal, the covariance matrix of the 

demodulation matrix is, in general, non-diagonal because matrix inversion 
.^ | is a nonlinear operation. This propagates through the demodulation process 

and induces correlations on the inferred Stokes parameters. © 2008 Optical 
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Introduction 



The majority of the information that we have obtained during the last years about the mag- 
netism of the Sun and other astrophysical objects is based on the analysis of the polarization 
of the light. The polarization state of a light beam is typically described, using the Stokes 
formalism, by the so-called Stokes vector: 

S=(I,Q,U,V) T , (1) 



where / refers to the total intensity of the beam, V describes its circular polarization prop- 
erties while Q and U are used to described linear polarization. When the light beam passes 
through optical elements or, in general, any medium that modifies its polarization properties, 
the emergent Stokes vector can be related to the input one by the following linear relation: 

S out = MS in . (2) 

The 4x4 matrix, M, is the so-called Mueller matrix and it can be used to unambiguously 
represent any optically passive medium. It is important to note that these matrices have to 
obey certain properties so that they have physical meaning [1,2]. 

The analysis of polarization of a given beam is usually carried out with the aid of mod- 
ulation schemes. This is a requisite in the short-wavelength (optical) domain because no 
detectors that present the required polarization sensitivity are still available. Furthermore, 
the majority of the detectors that are used in the optical are only sensitive to the total 
amount of light or, in other words, to the total intensity given by Stokes I. This is not the 
case in the long-wavelength (microwave) domain, where it is relatively easy to build detectors 
that analyze directly the polarization properties of the light beam. Modulation schemes have 
been built to overcome the difficulties when measuring the polarization state of light beams 
in the short- wavelength domain. Any modulation scheme consists of a train of optical devices 
that produce a known modification of the input beam so that the observed intensity in the 
detector is a linear combination of the elements of the input Stokes vector. Carrying out 
several of these measurements, each one with a different combination of optical devices, it 
is possible to infer the input Stokes vector using a demodulation procedure. Assuming that 
the Mueller matrix of the j-th combination of optical devices is given by Mj, the detected 
intensity in each case is given by: 

J° ut = [Mj-lorf + [MjUS? + [M^S? + [MAnS?, (3) 

where we have used the standard notation S = (So, Si, S2, S3) 71 for the Stokes vector, which 
will be used extensively in the rest of the paper. By putting several measurements together, 
the modulation scheme can be expressed as the following linear system: 

put = QS in , (4) 

where each row of the matrix O is built with the first rows of the different Mueller matrices of 
the different combinations of optical elements used in the modulation scheme. Therefore, the 
matrix O has dimensions iV x 4, with N the number of measurements. Note that, neither 
the vector I out is a Stokes vector nor the matrix O is a Mueller matrix. The vector I out 
consists of different intensity measurements, while the matrix O does not fulfill, in general, 
the conditions described by [1] and [2] to be considered a Mueller matrix. 



Several different modulation techniques have been developed. One of the most widespread 
methods is to use a temporal modulation with a quarter-wave plate and a linear polarizer, 
with the fast axis of the quarter-wave plate being set to preselected angles [3,4]. This approach 
has the advantage that the different polarizations states are observed on the same pixel. 
However, modulation has to be carried out very fast in order not to be limited by atmospheric 
fluctuations. Another possibility is to use spatial modulation in which the measurements 
are carried out simultaneously by using a beamsplitter instead of a linear polarizer [5,6]. 
The advantage is that atmospheric fluctuations are avoided but potential problems appear 
because the optical path and the detector response are not the same for all the polarization 
states. The most advanced polarimeters now use a combination of both approaches in the 
so-called spatio-temporal modulation [7, 8] . 

Starting from Eq. (T4J), the input Stokes vector can be obtained by solving the previous 
linear system of equations, so that: 

gin = D jout ( 5 ) 

where the matrix D is the inverse (or Moore-Penrose pseudoinverse in the general case that 
more than four measurements are carried out to infer the input Stokes vector) of the matrix 
O and has dimensions 4 x N. In order to measure the four Stokes parameters, the minimum 
number of measurements is 4. In such a case, the matrix D is unique and can be calculated 
as the standard inverse of O. However, it is also possible to carry out more measurements 
than Stokes parameters. In this case, the linear system of Eq. (j4j) is overdetermined and, 
in general, no solution that satisfies simultaneously all the equations exist. Recently, it has 
been shown that it is still possible to choose a solution if we seek for the one that maximizes 
the efficiency of the modulation-demodulation scheme [9]. Such efficiency for each Stokes 
parameter, represented with the vector e, is obtained from the demodulation matrix [9,10]: 

e;=U|XJ • (6) 

In this case, it can be shown that the optimal demodulation matrix is given by the Moore- 
Penrose pseudoinverse [11,12]: 

D=(0 T 0) _1 T (7) 



The pseudoinverse can also be found with the singular value decomposition of the modulation 
matrix. Note that, while DO = 1 in the general case, OD = 1 is only valid when N — 4 (1 
is the identity matrix). 

We present a theoretical calculation of how the error propagates in the demodulation 
process. This is an important issue that has been partially treated in the past, either con- 
sidering how the uncertainty in the measurement of I out propagates to the demodulated 



Stokes vector when the modulation matrix is perfectly known [13, 14] or taking into account 
uncertainties in the knowledge of the modulation matrix, though the full covariance matrix 
is not obtained [15, 16]. Many papers also deal with how error is propagated in non-ideal 
Mueller matrix polarimeters [17,18] while others consider the optimization of such polarime- 

ters [19-21]. 

2. Error propagation 

In practical situations, the modulation matrix has to be measured once the optics and the 
modulation scheme have been defined [14,16,17,20,21]. As a consequence, one expects the 
elements of the modulation matrix to have some uncertainty produced by the measurement 
procedure. Additionally, uncertainties can also be found when the rotation of the fast axis 
of the retarders or of the linear polarizers has some repeatability problems [22]. In this case, 
although the average modulation matrix can be known with great precision, random devi- 
ations from the average modulation matrix can occur. These random deviations induce a 
modification on the modulation scheme that directly affects the inferred polarization prop- 
erties of the input light beam. For these reasons, it is fundamental to characterize the error 
propagation through the demodulation process and how they affect the measured Stokes 
profiles. In general, since the matrix inversion (or pseudoinverse) is a nonlinear process, one 
should expect to find a full covariance matrix for the elements of the matrix D, with non-zero 
non-diagonal elements. These non-diagonal elements account for the statistical correlations 
that appear between elements of the demodulation matrix after the inversion process. More 
importantly, a full covariance matrix has to be expected even in the case that the elements 
of the modulation matrix are not correlated (i.e., cov(Oij,Oki) oc SikSji, with 8ik the Kro- 
necker delta). Note that a diagonal covariance matrix for the matrix O comes out when the 
measurements of the elements of the modulation matrix are statistically independent, some- 
thing that can be assumed in some cases if the measurements of the elements of the matrix 
are carried out with proper calibration. However, this is an issue that has to be carefully 
analyzed for each case. 

2. A. Analytical approach 

The full covariance matrix for the demodulation matrix can be obtained analytically following 
the standard error propagation formulae. The starting point is to consider that elements of 
the demodulation matrix are known nonlinear functions of the elements of the modulation 
matrix: 

D af3 = D aP {O vj ). (8) 



The most general error propagation formula for such operation reads [23]: 



cov(D a p, D ab ) = J2 -7^h-^7V- COY (°ij> °ki)- (9) 

ijki dUi i dUkl 



This covariance matrix gives information about the variances of the elements of the de- 
modulation matrix (diagonal elements), as well as covariances between those elements (non- 
diagonal elements). Note that, for a4xiV matrix, the covariance matrix is AN x 47V. If the 
covariance matrix of the modulation matrix is diagonal (i.e., cov(Oy, Oki) = cr(0)fj5ikSji, 
with 5ik the Kronecker delta), the previous equation simplifies to: 

„ v(D ^) = E^^(o)?, do) 

ij l j l j 

As already presented above, this turns out to be a good approximation if the elements of the 
modulation matrix are carefully measured. However, our approach is general and can cope 
with non-diagonal covariance matrices as will be shown below. 

The derivatives that appear in Eqs. (Q and (fTUl) can be calculated, in the most general 
case, starting from Eq. (|7j). However, we distinguish the cases of a square modulation matrix 
from the general case, because the analytical procedure can be largely simplified. In any 
case, we will show that the general expressions are equivalent to those of the particular case 
of a square matrix O. 

2.A.I. Square modulation matrix 

Although it is possible to calculate directly the derivatives from Eq. (JTj), it is more ad- 
vantageous to use the fact that the matrix product between O and D commute, so that 
DO = OD = 1. In this case, it is easy to show that [24]: 

dD aP 



dO l3 



D ai D 3P . (11) 



Therefore, the full covariance matrix of the matrix elements of the demodulation matrix can 
be calculated by substituting Eq. (flTT) into Eq. ([9]): 



cov(D a/3 , D ab ) = Y^ DuiDjpDakDibCovtOij, O kL ). (12) 

ijkl 

In the particular case of a diagonal covariance matrix for O, the previous equation simplifies 
to: 

cov(D af3 , D ab ) = £ D m D jP D ai D jb a{0)l. (13) 



2. A. 2. Non-square modulation matrix 



When N is larger than four, it is not possible to follow the previous approach because O 
and D do not commute. However, inserting the intermediate matrix A = O t O, using the 
fact that D = A~ x O T and that the matrix A is always invertible [9], the derivative can be 
expressed, after some algebra, as: 



9D aP ! ^ dA a l k 

- A aj dm + 2^ Opk 



dOij 



dO VJ • 



(14) 



The derivative of the elements of the A x matrix with respect to the elements of O can be 
calculated using the chain rule: 



dA~l 



dO, 



ak 



E 



dA-l dA r 



ak 



r) A f)f) ■ ■ 

ran KJ - n -mn wx ^ij 



The definition of the matrix A allows us to obtain the following derivative easily: 

oA mn 



80 



(15) 



(16) 



<:l 



The derivative of the elements of A x with respect to the elements of O can also be calculated 
easily from Eq. (ITT]) because A is a square matrix, that commutes with its inverse. Therefore: 



dOu 



A' 1 A' 1 

A ai A jp ■ 



(17) 



Substituting Eq. (ITTj) into Eq. (TT4l) and after some algebra, we end up with the following 
expression for the derivative: 

dD aP 



dO^ 



A a j S/3i — A a j 2^ DnfjOin — D ai Djp. 



This expression reduces to Eq. (ITTI) if O is square because the product of D and O commutes, 
so that: 

E D nf3 O in = Y^ O in D n/3 = 5ip, (19) 

n n 

After substitution in Eq. (ITBI) . we recover Eq. (TTTT) . 

The final expression for the covariance of the demodulation matrix is obtained by plugging 
Eq. (HE]) into Eq. © [or Eq. ([10])]: 



cov(D a/3 ,D ab ) 



J2cov(O i:j ,O kl ) 

ijkl 

A a j 5pi — A a j 2_^ D n pOi n — D ai Djp 

A-ak $bk — A ak 2_^ DmbOkm ~ Aifc Ab 

m 

6 



(20) 



When both the modulation matrix O and the measured intensities I out present uncertain- 
ties, the resulting uncertainty in the demodulated Stokes vector has contribution coming from 
both origins. Using standard error propagation formulae applied to Eq. (jSJ), the covariance 
matrix can be written as: 

arm Qjin 

cov(S?,S?) = V ' J cov(D a p,D ab ) 

+ E^^cov ( /r,/r). (2i) 

The first contribution takes into account the uncertainty in the knowledge of the modulation 
matrix while the second contribution takes into account the uncertainty in the measurement 
of the intensities arriving to the detector [13,25]. In the field of optimization of polarimeters, 
a suitable norm of this very last term (often neglecting the non-diagonal covariances in the 
vector I out that we include here) is the chosen one to measure the efficiency of a polarimeter 
[13,19,25,26]. Particularizing to our problem and calculating the partial derivatives, we end 
up with: 

cov(S*,£f) = ^CVcovfA,^) 

a/3 

+ J2D ik D jlC ov(ir,in- (22) 

ki 

The covariance matrix cov(I out ) is typically assumed to be diagonal, since no statistical 
dependence is assumed between consecutive intensity measurements in the detector, so that 
cov(/£ ut , I° nt ) oc Shi- Since this is not generally the case for the quantities cov(Di a ,Djp), a 
non-zero (in general) covariance between different Stokes profiles appears as a consequence 
of some degree of correlation in the demodulation process. 

In principle, it is possible to diagonalize the covariance matrix given by Eq. (122J1 to obtain 
its principal components. Such principal components define the directions along which the 
correlation between different Stokes parameters is minimized. Since the principal components 
define a new reference system, it is possible to rotate the original reference system in which 
the four-dimensional Stokes vector is defined. This rotation minimizes in a statistical sense 
the cross-talk. This could be of interest when synthetic calculations have to be compared 
with polarimetric observations. In this case, in order to minimize the effect of cross-talk, it 
could be advantageous to compare the observations with projections along the eigenvectors 
of the synthetic calculations. 

A quantity that is also affected by uncertainties in the knowledge of O is the efficiency 
defined in Eq. (jSJ). Applying the error propagation formula, we get: 

cov(e Q , e p ) = Y, 7^7^T cov ( 0i i' ° kl )- ( 23 ) 

ijkl dU v d(J ki 



If we assume that the covariance matrix of the elements of the modulation matrix is diagonal, 
we get: 

cov( ea , e ,) = £j^J^a(0)^ (24) 



The derivatives can be calculated using the definition of the efficiency and the chain rule: 

Z^c,n an..' W 



where 

<9e 



« 



90 y j^ dD pq dO i:j 



D aq 5 pa n^ 2 [J2 D 2 al ) . (26) 



d£ M v , 



Plugging this expression into Eq. ([25]) . we can simplify it to read: 



E^ E^T^- (27) 



The derivative dD aq /dOij is obtained from Eq. (ITS1) in the general case and from Eq. ( TTT1) 
in the case of a square modulation matrix. 

According to the previous results, when the covariance matrix of the modulation matrix 
is diagonal and the variance is the same for all the elements, all the covariance matrices 
calculated in this section are proportional to the variance of the elements of O. As a conse- 
quence, reducing an order of magnitude the uncertainty in the elements of the modulation 
matrix reduces an order of magnitude the uncertainty in the demodulation matrix and in 
the inferred Stokes vector. 

2.B. Monte Carlo approach 

The previous analytical approach assumes that the error propagates normally. This is typi- 
cally a good approximation when the modulation matrix is far from singular. In other words, 
the rows of the modulation matrix have to be as linearly independent as possible. However, 
since the inversion is a nonlinear process, it would be possible to obtain distributions in the 
demodulation matrix (and as a consequence, in the inferred Stokes parameters) that are far 
from Gaussians. In order to verify this issue and also to test that the previous analytical 
approach gives the correct answer, we have carried out a Monte Carlo experiment. For sim- 
plicity, we assume that the measurement uncertainty is the same for all the matrix elements 
of the modulation matrix and that it is characterized by the variance a 2 , so that the full 
variance matrix is given by: 

a(0)l = a 2 . (28) 

This selection makes the presentation of the results easier, but the previous approach is 
general and can cope with non-diagonal covariance matrices. The Monte Carlo experiment, 



given an initial modulation matrix, generates N instances of such a matrix with added 
uncertainties: 

6, = + 74 <t(0), (29) 

where 7* is a normally distributed random constant with zero mean and unit variance. For 
each modulation matrix, we obtain the optimal demodulation matrix following Eq. (J7J). 
Finally, the statistical properties of the demodulation matrix are characterized by the full 
covariance matrix, given by: 



cov(A;,£>fcz) = E[D y D fcI ] - E[D tj ]E[D 



ki\ 



(30) 



where E[x] stands for the expected value of the x random variable. In such a Monte Carlo 
approach, we estimate it using the sample mean, so that E[x] = (x). 

3. Illustrative examples 

The previous analytical approach is applied to several modulation matrices obtained from 
the literature and that represent different typical modulation schemes. We analyze the un- 
certainties in the demodulation matrix and in the inferred Stokes parameters using the error 
propagation formulation presented in the previous section. For illustration purposes, the 
results are also accompanied in some cases with plots obtained using the Monte Carlo ap- 
proach. One of the reasons for that is to demonstrate whether the resulting distributions 
can be correctly assumed to be Gaussians. The Monte Carlo results have been obtained 
using N = 60000. We assume that the variance that characterizes the uncertainty in the 
measurement of the modulation matrix is a 2 = 6.25 x 10~ 6 (the standard deviation is equal 
to a = 2.5 x 10~~ 3 ), which is a reasonable value for a standard calibration of the modulation. 

3. A. Square modulation matrix 

The first example is representative of a square modulation matrix in which four measure- 
ments are carried out to infer the four Stokes parameters. The matrix is that used by the 
Tenerife Infrared Polarimeter (TIP; [27]) as presented in [9]. Although this matrix is not the 
experimental one for TIP, its presents the desired structure and can be used to gain some 
insight on the properties of the error propagation. The chosen modulation matrix is: 



O 



TIP 



1 0.47 -0.68 0.48 

1 -0.91 -0.19 -0.13 

1 -0.11 0.57 0.72 

1 0.68 0.27 -0.58 



(31) 



D 



TIP 



0.19 


0.31 


0.18 


0.31 


0.31 


-0.66 


-0.03 


0.37 


-0.72 


-0.24 


0.67 


0.29 


0.36 


-0.34 


0.60 


-0.61 



(32) 



The optimal demodulation matrix, obtained from the application of Eq. ((Tj), is just the 
inverse of O: 

As shown in Fig.fU the distribution of the elements of the demodulation matrix are Gaussians 
whose standard deviation can be obtained from the diagonal elements of the matrix defined 
in Eq. Q3J: 

" 1.14 1.08 1.19 1.07 

1.81 1.71 1.88 1.70 

2.34 2.22 2.43 2.20 

2.19 2.07 2.27 2.06 



<x(D 



TIP, 



X 10' 



(33) 



Although we only present the diagonal elements of the covariance matrix, note that the 
covariance matrix is a full matrix in which all the elements are different from zero. 

As a direct consequence of the uncertainties in the modulation matrix, the efficiency of 
the proposed scheme presents uncertainties. They can be calculated using Eq. ( l24"l) . so that 
the average value and their standard deviations are: 



e TIp = (0.969,0.612,0.473,0.506) 
o-(e T ip) = (1.38, 1.26, 1.25, 1.27) x 10" 3 , 



(34) 



where non-diagonal elements in cov(e a , e@) are non-zero, but typically much smaller than the 
elements in the diagonal. 

Now we investigate the propagation of uncertainties to the inferred Stokes parameters. To 
this end, we choose an initial light beam defined by S m = (1, 10~ 3 , 10~ 3 , 10~ 3 ). Such Stokes 
vector is representative of what one would observe in relatively low magnetic flux regions of 
the solar surface. The application of Eqs. (1221) and ([TBI gives the following covariance matrix: 



covfS 1 



1.66 -0.21 0.01 -0.76 

-0.21 4.17 0.14 0.59 

0.01 0.14 6.99 0.32 

-0.76 0.59 0.32 6.09 



x 10' 



(35) 



This result has been obtained assuming that there is no uncertainty in the measurement of the 



intensity arriving to the detector, so that cov(J 



out rout"! 

k j 1 l ) 



0. In so doing, we isolate the effect 



of the modulation on the inferred Stokes vectors. The results of the Monte Carlo simulation 
are shown in Fig. [2], demonstrating that the values are quasi-normally distributed around 
the original value with a dispersion that is given in each panel of the plot. Note that the 
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standard deviations are large, almost of the order of the standard deviation of the precision in 

the measurement of the modulation matrix. Note also that the largest value of cov(5^ n , SJ 1 ) 

are on the diagonal of the matrix, but sizable non-diagonal elements also appear, indicating 

a certain degree of correlation between the inferred Stokes parameters. For instance, the 

value cov(S'o n , S™) = —0.76 x 10~ 6 indicates that when the inferred Stokes i" (S™) increases, 

the inferred Stokes V (S™) tends to systematically decrease. This can be understood as a 

cross-talk between all the Stokes parameters induced by the demodulation process due to the 

special structure of the modulation matrix. The two-dimensional distribution of the inferred 

Stokes I and V is shown in the left panel of Fig. [31 In specific cases, part of this cross-talk 

can be corrected for easily based on physical arguments. This happens, for instance, when 

observing the Stokes profiles induced by the Zeeman effect in magnetized regions of the solar 

surface with the aid of spectropolarimeters. In such a case, one can assume that Stokes Q, 

U and V are zero away from the spectral line and carry out a correction of the cross-talk 

from Stokes I to Stokes Q, U and V. 

The diagonalization of the cov(S m ) matrix gives the following eigenvector matrix ordered 

in row format: 

0.986 0.042 -0.012 0.160 

0.002 0.963 -0.018 -0.268 

-0.057 0.125 0.913 0.385 

-0.156 0.234 -0.408 0.868 



V 



(36) 



Each row represents the linear combination of Stokes parameters that one infers due to 
the induced cross-talk contamination. Note that the weight of one of the Stokes parameters 
is larger in each row, but the contamination from the other parameters is still large, a 
consequence of the chosen modulation matrix. 

3.B. Non-square modulation matrix 

We present here results for a typical non-square modulation matrix. We have chosen the one 
belonging to the Advanced Stokes Polarimeter (ASP; [28]) as presented in [9]. As it happens 
for the TIP matrix, this is probably not the exact matrix used in the polarimeter but is 
representative of what happens in a scheme in which more than four measurements are used 
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to obtain the four Stokes parameters. The modulation matrix we choose is: 



O 



ASP 



1 


0.77 


0.41 


-0.36 


1 


-0.06 


0.41 


-0.86 


1 


-0.06 


-0.41 


-0.86 


1 


0.77 


-0.41 


-0.36 


1 


0.77 


0.41 


0.36 


1 


-0.06 


0.41 


0.86 


1 


-0.06 


-0.41 


0.86 


1 


0.77 


-0.41 


0.36 



(37) 



After calculating the demodulation matrix and its covariance matrix, we can also calculate 
the diagonal of the covariance matrix for the efficiency which, transformed into standard 
deviations, give: 



e ASP = (0.760,0.415,0.410,0.659) 
<r(e A sp) = (1.37,0.94,0.88,0.88) x 10" 3 , 

and the full covariance matrix for the inferred Stokes parameters using S 1 
(1, 10~ 3 , 10~ 3 , 10~ 3 ) as the input Stokes vector: 



(3J 



covfS 1 



1.35 


-1.61 








-1.61 


4.54 














4.65 














1.80 



x 10" 



(39) 



Note that the elements of the diagonal of this covariance matrix are smaller than for the 
case of only four measurements (except for the case of cov(S'" 1 , Sf 1 )), probably induced by 
the larger number of measurements carried out. Furthermore, it is important to point out 
that this modulation scheme induces no correlations between S™ (Stokes U) and S l % (Stokes 
V) and any other Stokes parameter. On the contrary, there is a large correlation between 
Stokes / and Stokes Q, something that can be clearly seen in the right panel of Fig. [31 

The diagonalization of the previous matrix gives the following eigenvector matrix (in row 
order) : 

0.92 0.39 

-0.39 0.92 

10 

1 



V 



(40) 



The demodulation process has induced a cross-talk between Stokes / and Q which can be 
represented by a rotation of ~ 23°. 
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3. C. Ideal modulation matrix 

It is instructive to present results for the following ideal modulation matrix: 



O 



IDEAL 



1 


1 








1 


-1 








1 





1 





1 





-1 





1 








1 


1 











1 



(41) 



The efficiency in this case is the maximum that can be reached, while the diagonal of the 
covariance matrix for the efficiency given as standard deviations is: 



cideal = (1, v3, v^, \/3) 
<r(e IDE AL) = (1-02, 1.02, 1.02, 1.02) x 1(T 3 , 

and the full covariance matrix for the inferred Stokes parameters using S 1 
(1, 10" 3 , 10~ 3 , 10~ 3 ) as the input Stokes vector is: 



(42) 



covfS 1 



1.04 

3.13 

3.13 

3.13 



x 10" 



(43) 



The covariance matrix is diagonal with equal uncertainties in Stokes Q, U and V. These 
values are smaller than for the ASP example except for Stokes V, where the uncertainty 
for the ideal modulation matrix is slightly larger. More important is the fact that, since the 
covariance matrix is diagonal, no correlation is found between the inferred Stokes parameters, 
so that no residual cross-talk induced by the modulation is induced in the demodulation. 

4. Concluding remarks 

We have presented analytical expressions for the calculation of the propagation of errors in 
the demodulation process when the modulation matrix is not known with infinite precision. 
This can happen when the modulation system has not been calibrated with enough precision 
or when the repeatability of the modulation system induces uncertainties in the modulation 
scheme. The formulae that we have presented allows polarimeter designers to calculate the 
errors in the demodulation matrix, the efficiency of the modulation scheme and the inferred 
Stokes parameters. They are simple to calculate and require only the knowledge of the mod- 
ulation matrix, together with its covariance matrix. We have pointed out the fact that, in 
general, since matrix inversion (or Moore-Penrose pseudoinversion) is a nonlinear operation, 
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non-zero non-diagonal covariances have to be expected in the demodulation matrix even if 
such non-diagonal correlations are not present in the modulation matrix. This has the impor- 
tant consequence of generating spurious correlations (cross-talk) between the inferred Stokes 
parameters. We calculate the induced cross-talk by diagonalizing the covariance matrix. The 
matrix of eigenvectors represent the reference system in which the cross-talk is minimized. 
We have illustrated these points with three different modulation matrices representing three 
different ways of measuring the polarization state of light beams. Each method has its own 
advantages and disadvantages. It is up to the polarimeter designer to choose the modulation 
scheme depending on the desired precision. We hope that the formulae present in this pa- 
per are of interest for improving the quality of modulation based polarimeters. Routines in 
FORTRAN 90 and IDIJU for the calculation of the formulae presented in the paper can be 
obtained after an e-mail request to the authors of this paper. 
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Fig. 2. Distribution of the inferred Stokes parameters when no uncertainty is 
assumed in the measurement of the I out vector but only in the measurement 
of the TIP modulation matrix, whose uncertainty is a = 2.5 x 10 -3 . The 
input Stokes vector was S in = (1, 10~ 3 , 1CT 3 , 1CT 3 ). The distributions are close 
to Gaussians centered on the original values. The standard deviation of the 
inferred value is given in each panel. Note that the errors are roughly similar 
or slightly smaller than the uncertainty in the modulation matrix. Note also 
that correlations between the elements of the inferred vector are also expected 
and are produced by the non-diagonal covariance matrix of the demodulation 
matrix. 
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Fig. 3. Left panel: two-dimensional distribution obtained from the Monte Carlo 
simulation for S 1 ™ and S 1 ™ using TIP's modulation scheme, i.e., the inferred 
Stokes I and V, respectively. The plot shows the presence of a small correlation 
between the two results, a direct consequence of the appearance of correlations 
during the inversion process of the modulation matrix. Right panel: the same 
result but for the ASP case and for S™ and Sf 1 , i.e., the inferred Stokes / and 
Q. Note that the correlation between both is even larger, as indicated by the 
respective covariance matrices given by Eqs. (|35|) and (|39|) . 



